Design  Optimization 


Different  techniques  to  make 
designs  better! 
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What  is  an  Optimization 
Problem? 

“Optimization  is  a  process  of  selecting  or 
converging  onto  a  final  solution  amongst 
a  number  of  possible  options,  such  that  a 
certain  requirement  or  a  set  of 
requirements  is  best  satisfied.  ” 

l.e.,  you  want  a  design  in  which  some  quantifi¬ 
able  property  is  minimized  or  maximized 
(e.g.,  strength,  weight,  strength-to-weight  ratio) 
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Some  Types  of  Optimization 
Problems 


•  Continuous  real  parameter(s): 
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Some  Types  of  Optimization 
Problems 


•  Combinatorial: 


Geometric: 


Waste  ratio  =  17.5% 


Figures  from:  K.  Lee,  “Principles  of  CAD/CAM/CAE  Systems,”  Addison- Wesley,  1999 
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Some  Types  of  Optimization 
Problems 


•  Structural 


(b) 


E.g.,  Minimize  weight,  maximize  strength 


Figure  from:  K.  Lee,  “Principles  of  CAD/CAM/CAE  Systems,”  Addison-Wesley,  1999 
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Formulating  the  Problem 


To  have  a  set  of  possible  solutions,  the  design 
must  be  parameterized.  The  objective  function 
must  be  defined  in  terms  of  those  parameters. 

Formulation  steps: 

a)  Identify  design  variables 

(e.g.,  engine_displacement,  pistonjdiameter) 

b)  Define  objective  function 

(e.g.,  maximize  power_to_weight_ratio) 

c)  Identify  constraints 

(e.g.,  1  inch  < piston_diameter  <12  inch) 
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Formulating  the  Problem 
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Formulating  the  Problem 
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Formulating  the  Problem 


Mathematically,  you  need  to  select: 


a)  Design  Variables  (vector  X,  solution  X  e  R " ) 


V 

'5.4 

e.g.,X  =  < 

x2 

•,x*  =  . 

6.0 

*3. 

9.8 

This  means  X*  is 
a  vector  of  n  real 
numbers. 


b)  Objective  Function  (F(X)) 

X  e  Rn  so  that  F(X*)  =  min  F(X) 


In  other  words,  the  solution  is  the  vector  of  real  numbers  X*  for  which 
F(X)  is  minimum. 
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Formulating  the  Problem 


c)  Constraints 
-  bounds 


X  <x  <x 


5 

Xl 

10 

e-g.,- 

5 

<• 

x2 

-<- 

10 

5^ 

X3 

10 

-inequality  Gf(X*)>0  i=l,2,...,m 

(e.g.,*;  +  x;>0  G1(X*)  =  j£+x'2) 

-equality  HJ(X')  =  0  j  =  \,2,...,q 

(e.g.,x"2-x‘3=0  H,(X')  =  x2-x 3*) 
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Formulating  the  Problem 


Constraints  act  as  a  guide  for  the  optimization 
problem. 

Bounds:  Are  direct  limits  on  the  values  a 
parameter  can  take  (e.g.,  5  <  x1  <  10.) 

Inequality:  Are  expressions  that  limit  the 
values  parameters  can  take 

(e.g.,  x1  -  x2  —  5  >  0.) 

Equality:  These  reduce  one  design  variable 
for  each  equality  constraint. 

(e.g.,  x1  —  x3  —  5  =  0.) 
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Structural  Optimization 


Structural  optimization  is  a  specific  class  of 
optimization  problems  that  uses  an  FE  analysis  as 
part  of  the  objective  function  or  constraints. 


Structural  optimization  involves  three  elements: 

1 .  Automatic  modification  of  structure/FE  mesh 

2.  FE  analysis 

3.  Optimization  algorithm 


MEEM4403  Computer-Aided  Design  Methods 


EM^wTech 


6 


Structural  Optimization 


Three  ways  of  modifying  structures: 

-  Parameter  optimization 

•  Concept:  change  solid  model  parameters. 

•  Procedure: 

1 .  Create  initial  solid  model. 

2.  Create  initial  FE  model  from  solid  model. 

3.  Execute  FEA. 

4.  Evaluate  objective  function  and  constraints. 

5.  Stop  if  design  is  optimal;  otherwise: 

a.  Change  dimension  or  parameter  in  solid  model. 

b.  Re-execute  solid  construction. 

c.  Re-mesh  FE  model. 

d.  Return  to  Step  3. 
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Structural  Optimization 


Shape  Optimization 

•  Concept:  Move  part  boundaries. 

•  Procedure  (automated): 

1 .  Create  initial  FE  model;  or  create  surface/solid  model  and  mesh 

2.  Execute  FEA. 

3.  Evaluate  objective  function  and  constraints. 

4.  Stop  if  design  is  optimal;  otherwise: 

a.  Move  node  in  FE  model;  or  move  control  point  on  surface  and  remesh 

b.  Return  to  Step  2. 

Topology  Optimization 

•  Concept:  change  density  of  material  regions  to  form  shape  and 

topology. 

•  Procedure: 

1 .  Initialize  densities  in  FE  model. 

2.  Execute  FEA 

3.  Evaluate  objective  function  and  constraints. 

4.  Stop  if  design  is  optimal;  otherwise: 

a.  Correct  densities 

b.  Return  to  Step  2. 
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Shape  Optimization 


Optimal  truss  design:  Node  locations  and  cross-section 
properties  are  the  design  variables. 


Figure  from:  K.  Lee,  “Principles  of  CAD/CAM/CAE  Systems,”  Addison- Wesley,  1999 
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Shape  Optimization 


The  locations  of  the  FE  nodes  or  the  B-Spline 
control  points  are  the  design  variables. 


Figure  from:  K.  Lee,  “Principles  of  CAD/CAM/CAE  Systems,”  Addison- Wesley,  1999 
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Topology  Optimization 


In  topology  optim¬ 
ization,  the  design 
variables  are  the 
amounts  of  material 
in  each  cell. 

Material  is  only 
added  where  it  is 
needed  to  carry 
loads. 


Figure  from:  K.  Lee,  “Principles  of  CAD/CAM/CAE  Systems,”  Addison-Wesley,  1999 
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Topology  Optimization 


Besides  allowing  Rectangular  hole  Rank  2  Artificial  (43) 

for  size  and 
shape  changes, 
topology  optim- 

.  I.  20x20  mesh  20x20  mesh  20x20  mesh 

ization  allows 
voids  to  appear 
or  disappear. 

40  x  40  mesh  40  x  40  mesh  40  x  40  mesh 


80  x  80  mesh  80  x  80  mesh  80  x  80  mesh 

Figure  from:  K.  Lee,  “Principles  of  CAD/CAM/CAE  Systems,”  Addison-Wesley,  1999 
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Topological  Optimization 


Optimal  Smart  Structure  Design  at  MTU! 


Simultaneous  optimization  of: 

•  conventional  material  structure, 

•  embedded  sensors, 

•  embedded  actuators,  and 

•  controller  design. 


Design  Solution 


Design  Problem 


constant 
electric  field 


Choosing  a  Solution  Method 


Calculus  based 


Guided  random 
search  technique 


X 


Enumerative  technique 


Numerical 

Symbolic 

Simulated 

annealing 

Genetic 

algorithm 

1 

Gradient-based 


Gradient-  and 
Hessian-based 


E.g.,  Steepest 

E.g.,  Newton's 

descent  method 

method 

Figure  from:  K.  Lee,  “Principles  of  CAD/CAM/CAE  Systems,”  Addison- Wesley,  1999 
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Symbolic  Solution 


The  symbolic  solution  is  found  by  solving 

dF 


dX; 


—  0  such  that 


d2F 


d2F  d2F 


dxtdx2  J  dx2  dx2 


<  0  (for  n  =  2) 


•E.g.,  F(X)  =  xf  +  (x2-5f 

—  =  2x,=0,  —  =  2(i,-5)  =  0 

oxl  ox2 


x,  =  0,  x2=5 
(o)2  -  2  •  2  <  0 


•  However,  most  engineering  optimization  problems  have 
objective  functions  that  are  too  complicated,  with  many 
constraints,  and  with  too  many  design  variables  to  make 
this  method  tractable. 


MEEM4403  Computer-Aided  Design  Methods 


EM^mTech 


Gradient-Based  Solving 


•  Gradient-Based  methods  are  iterative.  They 
choose  a  better  solution  by  following  the  down¬ 
ward  slope  of  the  curve/surface  given  by:  <!!L 


X*  x 

•  The  problem  of  oscillations  in 
the  steepest  descent  method 
is  solved  by  Newton’s  method. 

Figure  from:  K.  Lee,  “Principles  of  CAD/CAM/CAE  Systems,”  Addison-Wesley,  1999 
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Local  Minima 


•  Gradient-based  methods  do  not  work  well 
when  there  are  several  local  minima: 


•  The  Simulated  Annealing  and  Genetic 
Algorithm  methods  were  introduced  to  solve  this 
problem. 

Figure  from:  K.  Lee,  “Principles  of  CAD/CAM/CAE  Systems,”  Addison- Wesley,  1999 
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Simulated  Annealing 


•  Concept:  when  temperature  is  high,  solution  will 
jump  around  a  lot  but  as  the  temperature  cools 
the  solution  will  settle  into  a  local  minimum 

•  Procedure: 

1 .  Set  initial  temperature  and  initial  guess  at  solution. 

2.  Evaluate  objective  and  constraints. 

3.  Stop  if  minimum  temperature  reached;  otherwise: 

a.  Take  a  random  jump  to  a  neighboring  solution. 

b.  Evaluate  objective  and  constraints. 

c.  If  solution  is  better  then  accept  it;  otherwise: 

a.  Generate  a  random  number 

b.  If  the  random  number  is  less  than  the  temperature 
then  accept  the  new  solution 

4.  Decrease  the  temperature 

5.  Go  to  step  3. 
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Simulated  Annealing 


•  As  the  “temperature”  cools,  jumps  to  higher 
values  of  F(x)  are  less  likely  to  be  accepted. 


F(x) 
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Genetic  Algorithm 


•  Concept:  have  a  large  set  of  solutions; 
combine  the  best  solutions  to  get  even  better 
solutions 

•  Procedure: 

1 .  Create  a  population  of  initial  random  solutions. 

2.  Evaluate  objective  and  constraints  for  each. 

3.  Create  new  solutions  for  next  generation  by: 

a.  Randomly  crossing  two  solutions 

b.  Randomly  creating  mutated  solutions 

4.  Evaluate  objective  and  constraints  for  each  new  solution 

5.  Get  rid  of  the  worst  solutions  so  that  the  total  number  of 
solutions  stays  the  same. 

6.  Increment  the  generation  number. 

7.  Stop  if  the  maximum  number  of  generations  is  reached. 
Otherwise  go  to  step  3. 
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Genetic  Algorithm 


•  In  genetic  algorithms,  the  parameters  in  X 
are  represented  in  a  binary  form: 

5  4.85  13 

1  I  ♦ 

Parameter  #1  Parameter  #2  Parameter  #3 

\  1  / 

101  I  0111100101  I  01101 

♦ 

101011110010101101 


Figure  from:  K.  Lee,  “Principles  of  CAD/CAM/CAE  Systems,”  Addison-Wesley,  1999 
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Genetic  Algorithm 


New  solutions  are  added  by  modifying 
existing  solutions  randomly  using: 

1 .  Cross-over: 


New 
children 
with  gene 
segments 
from  both 
parents 


111111111111111111111 

000000000000000000000 


111111111111111111111 

ooooooooooooooooboooo 


111111111111111 100000 
000000000000000011111 


Figure  from:  K.  Lee,  “Principles  of  CAD/CAM/CAE  Systems,”  Addison-Wesley,  1999 
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Genetic  Algorithm 


2.  Mutation: 


These  bits 
chosen  randomly 
to  be  changed. 


After  each  iteration,  only  the  solutions  that 
best  meet  the  objective  function  are  kept. 


Figure  from:  K.  Lee,  “Principles  of  CAD/CAM/CAE  Systems,”  Addison- Wesley,  1999 
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Constraint  Handling 


A  common  way  of  handling  constraints  is  to 
introduce  penalty  parameters  (e.g.  pk)  as 
multipliers  that  modify  the  objective  function. 


(Optimize  D(x,  pk) 
instead  of  F(x).) 


Cl  x 

Figure  from:  K.  Lee.  “Principles  of  CAD/CAM/CAE  Systems,”  Addison-Wesley,  1999 
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